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Abstract 

Monte Carlo simulation performed on a lattice system of biaxial molecules possess- 
ing D<ih symmetry and interacting with a second rank anisotropic dispersion potential 
yields three distinct macroscopic phases depending on the biaxiality of the constituent 
molecules. The phase diagram of such a system as a function of molecular biaxiality 
is greatly modified when a transverse dipole is considered to be associated with each 
molecule so that the symmetry is reduced to Civ ■ Our results indicate the splitting of 
the Landau point i.e. the point in the phase diagram where a direct transition from 
the isotropic phase to the biaxial nematic phase occurs, into a Landau line for a system 
of biaxial molecules with strong transverse dipoles. The width of the Landau line be- 
comes maximum for an optimal value of the relative dipolar strength. The presence of 
transverse dipoles leads to the stabilization of the thermotropic biaxial nematic phase 
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at higher temperature and for a range of values of molecular biaxiality. The structural 
properties in the uniaxial and biaxial phases are investigated by evaluating the first 
rank and second rank orientational correlation functions. The dipole induced long 
range order of the anti-ferroelectric structure in the biaxial nematic phase, is revealed. 

1 Introduction 

In the past few years great interest has been paid to thermotropic biaxial nematic liquid 
crystals, whose existence was first predicted by Freiser pQ in 1970 from mean field theory. 
Shortly afterwards, a number of theoretical [2j [3] and computer simulation [U El [6] studies 
have been performed and phase diagrams consisting of three distinct macroscopic phases 
namely, the isotropic phase (J), the uniaxial nematic phase (Ny) (which may be, keep- 
ing the symmetry of the phase unchanged, of two types - prolate (N v +) or oblate (Ny-) 
nematic phases depending on the rod-like or disk-like nature of the molecules) and biax- 
ial nematic phase (Nb) were generated. On the experimental front, scientists have been 
engaged in search for truly biaxial nematic phase in thermotropic liquid crystals since its 
theoretical prediction. Although there have been a number of reports [7J [8] of experimental 
identification of thermotropic Nb phase since 1986, none of these claims proved to be cor- 
rect [11] . Recently there have been claims of observing the elusive thermotropic Nb phase 
for V-shaped molecules [121 US] and for tetrapode molecules [HI [T5]. These new findings 
have fuelled experimental and theoretical research [16] on biaxiality in thermotropic liquid 
crystals. 

The molecules forming liquid crystals always deviate from their assumed cylindrical 
symmetry and their usual structure can schematically be described as board-like (see Figure 
[T|). Usually this type of molecules give rise to the Nu phase as a consequence of the long 
range orientational order of the molecular long symmetry axis (w) while the short symmetry 
axes (u, v) remain uncorrelated except at short range. In this case the molecular long axes 
tend to be parallel, on average, along a single macroscopic direction called the principal 
director (n). At some state point long range orientational ordering is established for the 
molecular long axes as well as for the molecular short axes which results in a biaxial nematic 
phase and the corresponding average preferential alignments of the short axes occur along 
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Figure 1: Schematic diagram of a biaxial molecule. Three molecular symmetry axes are u, 
v and w. 

the two secondary directors (/, m). 

The ability of a material to form the Nb phase depends on the molecular biaxiality 
of the constituent molecules. In the case of isotropically averaged dispersive interactions 
(which we have used in the present work) there is a single molecular biaxiality parameter A 
[9J [10] which is related to the anisotropy of the polarizability tensor of the molecule. From 
the study of molecular field theory [TTJ it has been observed that for obtaining a stable 
Nb phase, A must have values within a very narrow range around the optimal biaxiality 
(Ac). The theoretical value of Ac is l/v^6 (~ 0.40825), which corresponds to the transition 
point from a system of prolate to oblate molecules [HI EU ELT] • For this value a second-order 
transition occurs directly between biaxial and isotropic phases and the corresponding point 
in the phase diagram is known as the Landau point. At this special critical point where 
two second order critical lines meet a first order phase boundary, Nb phases are expected 
to appear at the highest temperature and therefore, on the basis of theoretical analysis 
one may predict that to be a good candidate towards forming the Nb phase the molecule 
should have the biaxiality very close to Ac*. But it is almost an impossible task for the 
experimentalists to design mesogenic molecules with the desired biaxiality due to molecular 
flexibility. However, influence of other factors, such as polydispersity [IB], large transverse 
electric dipoles associated with the mesogenic molecules etc., enhances the stability of the 
Nb phase. The present scenario in context of the transverse dipole moments is elaborated 
below. 

For symmetric bent-core molecules (V-shaped) the optimal biaxiality factor Ac* = 1/ \/6 

3 



occurs for a single angle of 109.47° between the arms which is the tetrahedral angle. On 
the other hand Madsen et. al. [12] while investigating a class of bent-core molecules with 
different angles between the arms claimed to have observed a biaxial phase at an inter-arm 
angle as large as 140°. At this angle the molecular biaxiality factor is such that a Ny 
phase is expected. Although the work of Madsen et. al. [12] does not give a conclusive 
evidence of the occurrence of a biaxial nematic, it may be noted that the type of molecules 
they considered had a strong transverse electric dipole moment. To find if this could in 
some way make the job of finding the biaxial nematic phase easier, Bates [Tj5] performed a 
Monte Carlo simulation using a lattice model of bent-core molecules with transverse dipoles 
attached to the centre. In this work each molecule was assumed to have two symmetric arms 
and the included angle was varied. Each arm of a molecule interacts with each of the twelve 
arms of the six nearest neighbour molecules confined to a simple cubic lattice via the usual 
Lebwohl-Lasher potential [20]. In addition, the dipole-dipole interaction, which actually is 
of long range nature, was also confined to the nearest neighbours for simplicity and was 
taken to be of the — cos (3 type where (3 is the angle between the interacting dipoles. The 
phase diagram that Bates obtained was quite different from the usual one obtained while 
treating non-polar bend-core molecules where one gets a single Landau point for a bend 
angle of 109.47°. Instead, he observed that for large dipole moments the single Landau 
point gets transformed into a series of Landau points usually called a Landau line which 
covers a significant range of angles (and hence a range of molecular biaxiality) across which 
a direct Nb -h- / first order transition occurs. If this is really the case then the job of the 
experimentalists becomes a little easier. 

On the theoretical front there is a very recent work on bent-core molecules by Grzy- 
bowski and Longa [21] who used density functional theory (DFT) to analyze two- and 
three-segment Gay-Berne (GB) molecules with a single transverse dipole moment. They 
found that nonpolar models with two uniaxial arms has a single Landau point and the 
inclusion of the dipoles shifts the Landau point towards lower apex angles for the system 
containing two-arm molecules. For the three-segment molecules with large dipole moments 
a Landau line occurs and they observe that for an optimal dipole strength the Nb phase 
is most probable. In all cases these authors show that further strengthening of the dipolar 
interaction results in shrinking the distance between the isotropic phase and the biaxial 
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nematic phase which indicates that the intermediate uniaxial nematic phase becomes less 
stable. The study reveals that in the most likely situation for detecting the Nb phase each 
molecule should consist of three segments with a transverse dipole moment of about 3D. 

In this communication we present the results of a Monte Carlo simulation in a lattice 
model where biaxial molecules of D 2 h symmetry interact with nearest neighbours via an 
anisotropic dispersion potential. To study the effect of dipolar interaction on the phase 
behaviour of the system we have associated a transverse dipole moment with each molecule 
which reduces the symmetry of the molecules to C^. The dispersion interaction has been 
extensively used in mean field and MC studies [U [5j [17] to predict the phase behaviour 
of nematic liquid crystals. The case of symmetric V-shaped molecules considered in the 
work of Bates [T9] can be deduced from the dispersion potential we have used, as has been 
demonstrated in the work of Romano [22] and Bates and Luckhurst [16J. Although we have 
used the isotropically averaged out (over the intermolecular vector) form of the dispersion 
potential (which is short range in nature), the dipolar interactions have been considered 
in the full form and keeping in mind the long range nature of this interaction we have 
used the reaction field (RF) method [23j [24] in our calculation to improve the reliability 
of the results. Although the RF method is less rigorous than the Ewald summation (ES) 
method, it allows much computational advantages than the latter and the difference of the 
averages of thermodynamic quantities obtained by the RF and ES techniques are within 
the statistical uncertainty of the simulation [251 EHl 127J- We have compared our findings 
with the results of Bates [19] and Grzybowski and Longa [21] and these seem to be in good 
qualitative agreement. It may be noted that the dispersion potential we have used is a little 
more general than that in [19] since the former also includes the case of V-shaped molecules 
with non-symmetric arms. We have also investigated the details of the phase structures 
of this model by evaluating structural quantities like first-rank and second-rank correlation 
functions. 



2 The Model 

In our study we consider identical biaxial molecules with a transverse dipole, whose centers 
of mass are associated with a simple-cubic lattice. The interaction energy of two molecules 
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i and j is given by 

c,j r;f" • cj]' (i) 

where llfj sp is the pair potential obtained from London dispersion model [28| [29] and U^ 1 
is the dipole-dipole interaction. The dispersion term decays much faster (as 1/r 6 ) than 
the dipolar term (which varies as 1/r 3 ). We therefore, restrict the dispersion interaction 
to the six nearest neighbours. The orientationally anisotropic dispersion interaction explic- 
itly depends both on the mutual orientation of the two interacting molecules, and on their 
orientations with respect to the intermolecular unit vector (r^). By isotropically averag- 
ing over the intermolecular unit vector, the dispersion potential between two identical 
neighbouring molecules becomes 

u*» = + 2A[i4(%) + Ri (^,)] + 4A 2 i4(%)}. (2) 

Here = {<pij, 6ij, ipij} denotes the triplet of Euler angles defining the relative orientation 
of i th and j th molecules; we have used the convention used by Rose [5U] in defining the Euler 
angles, is the strength parameter which is assumed to be a positive constant (e) when 
the particles % and j are nearest neighbours and zero otherwise. R^ n are combinations of 
symmetry-adapted (D 2 h) Wigner functions 

R 2 00 = l^e- 1 - (3) 

Rl 2 = -^-sin 2 0cos2V> (4) 
/fi 

R 2 20 = 1— S i n 2 cos 2(f) (5) 

I 1 

i?22 = ^(1 + cos 2 6) cos 1(f) cos 2,tf) — - cos 9 sin left sin lip. (6) 

The parameter A is a measure of the molecular biaxiality and depends on the molecular 
properties. For the dispersion interactions, A can be expressed in terms of the eigenvalues 
(Pi) P2-, Pz) of the polarizability tensor p of the biaxial molecule 

, /3 P2-P1 

x -iii P3 - P2 - Pl - (7) 

The condition for the maximum biaxiality is P3 — p 2 = P2 — Pi > and A = Ac = 
(this self-dual geometry corresponds to the Landau point in the phase diagram). A < 
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Ac corresponds to the case of prolate molecules whereas A > A^ corresponds to oblate 
molecules. The V-shaped lattice model used by Bates [19] is equivalent to the present 
dispersion model and it can be shown [16] that both the models produce identical results 
when the temperature is rescaled properly by using 

AT* 

T* = VSM (R) 

(l-3cos 7 ) 2 1 ' 

and the molecular biaxiality parameter is given by 

[3 l + cos7 

A= V21-3cos 7 (9) 

where 7 represents the angle between the two symmetric arms of the V-shaped molecule. 
T* and Ty SM are respectively the dimensionless temperatures used in this model and in 
the V-shaped symmetric model studied by Bates. In general, the biaxiality parameter A 
depends both on the inter-arm angle and the relative mesogenic anisotropy of the arms for 
the symmetric as well as non-symmetric V-shaped molecules [16] . 

The fully anisotropic biaxial dispersion model has been studied for both two-dimensional 
and three-dimensional lattices by mean-field and MC methods [311 [22] . The MC simulations 
performed on a three-dimensional simple cubic lattice system of biaxial molecules interacting 
via the fully anisotropic dispersion model in one case and its isotropically averaged form 
in the other case suggest that both the models produce almost similar results [22]. We, 
therefore, opted the simpler form of the biaxial dispersion pair potential (Eq. ([2])) which 
has been investigated in previous studies and phase diagrams for a fairly wide range of 
molecular biaxiality have been generated [5j [T7J . 

The dispersion potential (Eq. (|2])) can also be expressed in Cartesian form [17] as 

U% sp = -e{[~(™ rW ^ 2 -^]-AV6^^ 

(10) 

where u, v and w are the three mutually orthogonal unit vectors along the molecular 
symmetry axes with the convention [32] that u and w are along the shortest and longest 
axes respectively (Figured]). Here the orientation of each biaxial molecule, in terms of the 
Euler angles, with respect to the laboratory frame {e x ,e y ,e z } is given by, 

u = (cos 9 cos cos ip — sin0sin'?/')e a . + (cos 9 sin cos ip + cos sin -^e^ — sin 9 cos ipe z , 
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v — — (cos 9 cos (f) sin ip + sin <p cos ip)e x — (cos 9 sin sin ■?/> — cos 4> cos t/>) e y + sin 6 1 sin ipe z , 
and 

iu = sin 9 cos + sin 9 sin 0e y + cos 9e z . 
The dipole-dipole interaction for a pair of point dipoles is given by 

Uff = 4^tK«* • «i) - 3(«< • ^)(«i " (11) 
Here is a unit vector along the inter-molecular vector connecting the centres of mass 
of molecules i and j of the dimensionless separation r^, /i is the common value of the 
transverse dipole moment, a is the nearest neighbour separation and Ui is the unit vector 
along the shortest molecular symmetry axis of the i th molecule. The long-range dipolar 
potential contributions to the total pair potential have been evaluated using the reaction 
field (RF) method. In this method each dipole is assumed to interact with all others within 
a spherical cavity of radius r c and beyond this cutoff radius the dipoles are considered to 
act as a dielectric continuum of dielectric constant e s producing a reaction field within the 
cavity. Therefore the dipolar pair potential within RF geometry can be written as 

tjw = jjdd _ 2(e, - l)/x 2 (Uj ■ Uj) 
ij 13 2e s + l 47re a 3 r c 3 ' 1 1 

for rij < r c and = for ry > r c . The truncation must be done for r c < L/2, L being 

the lattice dimension used in the simulation. The dielectric constant e s of the surrounding 

medium is usually chosen to be one of the two extremes, 1 for vacuum or oo for conductor. 

But a self consistent method [25] for estimating e s during the simulations is most appropriate 

in the RF method. We have used a parameter 

e * = T"^ ( 13 ) 
47re ecr 3 

which determines the relative strength of the dipolar interaction to the dispersion interac- 
tion. For a given value of e* the relative contribution of the dipolar energy in the total energy 
at the lowest temperature can be estimated by evaluating the ratio 2.676e*/(4.0 + 2.676e*) 
where the numerator is related to the negative of the ground-state energy (per particle) of 
a simple cubic dipolar lattice system [33J whereas the first term in the denominator is the 
negative of the lowest value (per particle) of the dispersion energy for the optimal value 
of the molecular biaxiality A = 1/ a/6. The dimensionless temperature used is defined as 
T* = k B T/e. 
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3 Computational Aspects 



To explore phase diagrams we performed a series of Monte Carlo (MC) simulations using 
the conventional Metropolis algorithm on a periodically repeated simple cubic lattice model, 
consisting of N — 40 3 particles. Simulations were run in cascade, in order of increasing 
dimensionless temperature T*. Equilibration runs were typically of 20 000 sweeps where 
each sweep consisted of iV MC steps (or moves) and were followed by a production run of 
20 000 — 30 000 sweeps (longer runs were used close to the transitions). 

In our case a Monte Carlo move was attempted by selecting a site at random and then by 
choosing one of the laboratory axes at random the molecule at that site was rotated about 
the chosen laboratory axis following the Barker- Watts method [33]. The dipolar energy of 
the molecules has been computed using the RF method, with cut off r c = 8 (in units of 
lattice spacing). The effect of the cut off radius on the simulation results have been studied 
in [351 [36]. I n our case further increase in r c produces almost indistinguishable results but 
increases the computation time significantly. For example if we choose r c = 10 the averaged 
quantities like order parameters and internal energy change by less than 0.1% compared to 
the results obtained for r c = 8. However, the CPU time was 62 hours for the former case 
and is about 28 hours for r c = 8 for simulation at a given temperature with 20 000 Monte 
Carlo steps per site on Intel Core i7 960 processors clocked at 3.2 GHz. The dielectric 
constant e s has been evaluated using a self consistent method [25]. In order to analyze the 
orientational order we have calculated the order parameters (R^n) following the procedure 
described in [37]. According to this, a Q tensor is defined for the molecular axes associated 
with a reference molecule. For an arbitrary unit vector w, the elements of the Q tensor are 
defined as 

Q afi (w) = ((3w a Wf}-6 a f ) )/2) (14) 

where the average is taken over configurations and the subscripts a and (3 label Cartesian 
components of w w.r.t the laboratory frame. The three Q tensors associated with three 
molecular symmetry axes are diagonalized (once after each MC sweep) and the eigenvalues 
and the eigenvectors obtained are then recombined to give the order parameters (R^n) w - r -t 
the director frame [37] . 

Of the four second rank order parameters only (Rq Q ) and (R% 2 ) are °f importance in 
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practice. (Rq ) is the usual uniaxial order parameter (P 2 ) which measures the alignment of 
the longest molecular symmetry axis with the primary director (n) and (R 22 ) is the measure 
of the alignment of the transverse molecular symmetry axes about the secondary directors 
(I, m). (-Rqo) is zero m an I phase and it is non-zero in the Njj phase as well as in the 
N B phase with maximum value 1 in the perfectly ordered state. On the other hand (-R| 2 ) 
is zero in the I phase and also in the Njj phase and is non-zero only in the Nb phase and 
increases towards its maximum value of 0.5 in perfectly ordered biaxial phase. 

We have also calculated the reduced specific heat per particle from fluctuations in the 
energy 

° v = Nk B T* (15) 
where E is the total energy of the system. 

Further details on the phase structure of the N B phase in presence of transverse dipoles 
can be found by analyzing structural quantities, like the first-rank orientational correlation 
function defined by 

0i(r) = <Pi(t*i • u 2 )) r , (16) 
and the second-rank orientational correlation function, 

g 2 {r) = {P 2 {u 1 ■ u 2 )) r . (17) 

The first, which is essentially just the average cosine of the angle between the transverse 
symmetry axes of molecule 1 and molecule 2 separated by a distance r, is used to determine 
the anti-ferroelectric structure of the biaxial phase, while g 2 (j') is useful in assessing the or- 
dering of the transverse molecular symmetry axes irrespective of their parallel or antiparallel 
arrangements. 



4 Results and discussion 

We present phase diagrams for five different values of the relative strength factor e* = 
0,0.4,0.8,1.2, and 1.6. The different phases are identified from the non-vanishing values 
of the relevant order parameters while the transition temperatures are obtained from the 
positions of the peaks in the temperature dependence of the specific heat curves (Cy) as 
well as from the order parameter curves and these were found to match with each other 
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Figure 2: The phase diagram as a function of dimensionless temperature T* and molecular 
biaxiality A (ranging from 0.15 to 0.6 in step of 0.05 along with the self-dual value of 0.40825) 
for a system of biaxial molecules interacting via second rank anisotropic dispersion potential 
only. The solid line represents the first-order transitions, while the dashed line represents 
the second-order transitions. 

within the errors present in the simulation. In the neighbourhood of a transition we have 
taken numerical data at temperature intervals of 0.025 and elsewhere the intervals were 
greater. To make a conclusion about the order of the phase transitions we have used the 
order parameter curves - for continuous vanishing of an order parameter the transition was 
taken as second order whereas a jump to zero or near zero value (which is a finite size effect) 
from a finite value was considered to be first order transition. These observations were in 
general supported by our results on specific heat plots. At a second order phase transition 
the peak heights are relatively smaller than those in a first order transition where the peaks 
are narrow and sharp. It may be noted that this method of identifying the phase transition 
temperatures and the order of the transition is only qualitatively correct. To ascertain these 
accurately one needs to perform more extensive MC simulation using, for instance, the finite 
size scaling methods which involve collecting data on systems of different size. We did not 
go into this exercise and therefore expect our results only to be qualitatively correct. 

In Figure |2] we have shown the phase diagram for e* = which means that only the 
isotropically averaged dispersion interaction (confined to nearest neighbours) given by Eq. 
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Figure 3: The phase diagrams of biaxial molecules interacting via the resultant potential 
with four different dipolar strengths e* = 0.4, 0.8, 1.2, 1.6. 



(J2J) has been considered. The result is same as obtained by Biscarini et al [5] and shows the 
J, Njj+, Nij- and the Nb phases. The biaxiality factor A has been varied from 0.15 to 0.6. 
The Nb <B- Njj transition is found to be second-order while Njj <-> I transition is first-order 
in nature and these are well known from previous studies [21 El El ED]- At the Landau 
point (i.e. for A = l/y/E ~ 0.40825) a direct Nb -H- / transition takes place and this is 
second-order in nature. In Figure El we have shown the phase diagrams when the long range 
interactions among the transverse dipoles have been considered in addition to the dispersion 
potential. The diagrams have been plotted for four values of the relative dipolar strength 
e* = 0.4,0.8, 1.2, and 1.6. When the strength of the dipolar interaction is weak (e* = 0.4 
and 0.8) the higher transition temperature j remains almost unaltered but the presence 
of (transverse) dipolar interaction stabilizes the Nb phase by raising the temperature T^ bNu 
to higher value. By increasing the relative dipolar strength, e* to 1.2 an indication of the 
formation of Landau line is observed and for strong dipolar interaction (e* = 1.6) the Nb ^h- 
/ transition occurs for a range of values of the biaxiality factor A (from 0.3 to 0.45) instead 
of a single value. Thus the Landau point (for e* = 0) turns into a Landau line in presence 
of strong transverse dipole moments (Figure E£0- Besides performing detail study for e* 
upto 1.6 we have performed a search for the extension of the Landau line for e* = 2.0,2.4 
and 2.8, where simulations were performed only to locate the two ends of the Landau lines. 
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We have observed, that the formation of Landau line starts at e* = 1.2 and the width of 
the Landau line i.e. range of the values of A showing the direct Nb -h- / transition, first 
increases and becomes maximum extending from A = 0.15 to 0.5 for e* = 2.0 at which the 
dipolar interaction is approximately 57% of the total energy at very low temperature. For 
further increase in the dipolar strength the width of the direct Nb -h- / transition becomes 
shorter and gets shifted towards higher values of A. For example when e* = 2.4, the Landau 
region covers the range of biaxiality parameter from A = 0.35 to 0.55 and for e* = 2.8 it 
extends from A = 0.55 to 0.7. 

We make a qualitative comparison of our findings with those of Grzybowski and Longa 
[21]. These authors used GB interaction between two-segment and three-segment molecules 
and also considered the effect of attaching a transverse dipole to each molecule. For two- 
segment molecules they do not find any Landau line resulting from the presence of the 
dipoles. This is in disagreement with our observation and with those of Bates [T9j . Moreover, 
we observe that after the formation of Landau line, with increase in the strength of the 
dipoles the extension of the Landau line first increases, reaches a maximum and then shrinks, 
always moving towards higher values of A. Grzybowski and Longa [21] for three-segment 
GB molecules find the same features in the extension of the Landau line. But initially this 
happens for lower values of the apex angle 7 (i.e. higher A) and beyond a certain value of 
the dipole strength the Landau line, which is of smaller extension, occurs at higher values 
of 7. Moreover, we observe that the most favourable situation (in terms of the width of the 
Landau line) for the occurrence of the Nb phase occurs for e* = 2.0. If one considers the 
lattice spacings (er) to be of the order of 5 A and uses the relation in Eq. ffT31) [38], then for 
a transition temperature of 473 K (which is close to that of the V-shaped molecule used in 
[i~2]). the estimated value of the dipole moment turns out to be 3.3 D which is close to the 
value of n obtained in [2T]. Here we have used the reduced transition temperature to be 1.5 
to estimate e. 

We now make a comparison of our results with those of Bates [19]. Bates observed 
that as the dipolar strength is increased, there is an increase in the Nb -h- Njj transition 
temperature. For e* = 0.5 the temperature at which the Nb phase is observed is essentially 
independent of 7 (or, A) and this is just below the temperature at which the Landau point 
is observed for e* = 0. This feature is similar to our phase diagram shown in Figure [3fc for 
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e* = 1.2. We are also inclined to state that for this value of the dipolar strength we find that 
the Landau line just begins to appear. Moreover, Bates' work shows that for e* = 1.0 there 
is a Landau line extending from 7 = 107° to 122° (the corresponding biaxiality parameter A 
being from 0.46 and 0.22) across which there is a first order Nb -h- I transition. This may 
be compared with our observation for e* = 1.6 where we find a Landau line extending from 
A = 0.30 to 0.45 across which there is a first order Nb -h- I transition. We however point 
out that a direct comparison of our results with those of Bates is not meaningful because 
of different features appearing in the phase diagrams for different values of the dipolar 
strengths. This is not surprising since the treatment of the long range dipolar interactions 
using the RF method rather than truncating it at the nearest neighbours is likely to make 
this difference. To state this more clearly, the qualitative features observed by Bates now 
appear at an enhanced value of e*. Also, while comparing the two sets of results, we need 
to rescale the temperatures using Eq. (jSJ). 

In Figures H] and [5] we show the plots of C v for A = 0.3, and A = 0.40825 respectively for 
the case of pure dispersion interaction (e* = 0) and for two other cases with relative dipolar 
strengths e* = 0.8 and 1.6. In Figure UK, which is for e* = and therefore corresponds to 
the phase point at A = 0.3 in Figure [21 we find that C v exhibits two peaks. The smaller 
peak at low temperature (T^ gN = 0.40) is for a Nb <-> N v transition and the sharper 
peak at higher temperature iT^ !uI = 1.10) is for the Nu <H- / transition. In Figure Hb 
(e* = 0.8) the two peaks of C v become closer due to the presence of dipolar interaction 
{T^ bNu = 0.85, T^ vI = 1.175). In FigureHt (e* = 1.6) a single sharp peak in C v is obtained 
which has a greater height and also occurs at a higher temperature {T^ gI = 1.30) indicating 
the formation of the Landau line. In Figure [5^ we find a single broad hump in the specific 
heat plot for A = 0.40825 and e* = which is similar to the earlier findings [5]. Figures 
[5b and [5fc exhibit the influence of the dipolar strengths on C v vs T* curves at the Landau 
point. We see that with the increase in dipolar strength the peak height increases and the 
transition temperature T^ gI shifts towards higher value. 

The temperature dependence of the second rank order parameters (Rqq) and (i?| 2 ) f° r 
e* = 0, 0.8 and 1.6 and for two different values of the biaxiality factor A = 0.3 and 0.40825 
are shown in Figure [6] and [7] respectively. In Figure [6^, which is for e* = we find that the 
uniaxial and biaxial order parameters (Rq ) and (-R22) represent order-disorder transitions 
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Figure 4: Specific heat per particle vs dimensionless temperature T* for the molecular 
biaxiality A = 0.3 and for three different dipolar strengths e* = 0, 0.8, 1.6. 
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Figure 5: Specific heat per particle vs dimensionless temperature T* for the molecular 
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at two well separated temperatures. The biaxial order parameter (P| 2 ) decays continuously 
almost from its maximum value of 0.5 to at the Nb <-> Njj transition (T^ bNu = 0.40) 
whereas the uniaxial order parameter (-Rqo) decreases with the increase in T* and jumps 
to zero at the Nu -h- / transition (T^ j = 1.10) which is an well known weakly first order 
transition. In Figure Eb (e* = 0.8) both the N B -H- Nu and Nu <-> / transitions get shifted 
towards higher T* and at the same time become closer due to the presence of dipolar 
interaction (T^ bNu = 0.85, T^ aI = 1.175). In Figure [6b (e* = 1.6) both the transitions 
transform into a single direct Nb <H- 1 transition at a higher temperature (T^ bI = 1.30). The 
behaviour of (Rq Q ) and (R% 2 ) with T* in Figure [7^ for e* = and A = 0.40825 qualitatively 
reveal the continuous nature of the direct Nb -h- / transition {T^ bI = 1.10). The increase 
of strength of the dipolar interaction to e* = 1.6 (Figure [7b) introduces a distinct change in 
character of the direct transition. In these figures sharp jumps in the values of (Rq Q ) and 
(R22) from a considerable finite value to nearly zero values are observed as T* is gradually 
increased (T^ z = 1.375). 

We have also calculated the first rank polar order parameter (P") = (tij ■ /) which is a 
measure of the degree of alignment of the transverse molecular symmetry axis u relative 
to the / direction of the director frame. We observe that this remains almost zero (within 
statistical error) over the entire temperature range for different values of the molecular 
biaxiality parameters and for all dipole strengths. Therefore, in any stable phase the system 
does not develop any spontaneous polarization. The ground state of the simple cubic system 
in presence of the dipolar interaction (Udd) alone is anti-ferroelectric and hence does not 
produce any net overall polarization. This is in agreement with the prediction made by 
Bates [19] for a system of real dipolar bent core molecules although with the Heisenberg form 
of dipolar interaction used by Bates is likely to yield a polar phase at low temperature for 
strong dipole moments. It may be pointed out that the dipolar interaction used in reference 
[T9"j is obtained by averaging out the anisotropic dipolar interaction over all orientations of 
the intermolecular vector for the nearest neighbours in a simple cubic lattice whereas our 
finding of the absence of a polar phase is obtained with the full form of the dipole-dipole 
interaction including its long range character. 

The orientational correlations between the transverse molecular symmetry axes are rep- 
resented by plotting the first-rank and the second-rank orientational correlation functions 
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gi(r) and g2(r), denned in Section [31 as a function of molecular separation r in Figure [8] 
for the molecular biaxiality A = 0.3. In case of pure dispersion interaction (e* = 0) the Nb 
phase at lower temperature (T* = 0.3) has a long-range behaviour of g2{f) (Figured). The 
long-range limit of g^O") (0.537) agrees with the square of the corresponding order parameter 
(which is 0.728), according to the relation 

\\mg 2 {r) = {(hu i .lf-\)Y. (18) 

r— >oo Z Z 

We see that gi(r) is zero for all values of r which shows, as can be anticipated, that the 
transverse molecular symmetry axis u has head tail symmetry in absence of transverse 
dipoles. 

At higher temperature (T* = 0.8), the long range order of the anti-ferroelectric structure 
in the Nb phase, induced by dipolar interaction (e* = 0.8), is observed from the pronounced 
oscillations of gi(r) as shown in Figure [8b. The positive value of g\{r) indicates a paral- 
lel arrangement of the pair of dipoles whereas its negative value shows their antiparallel 
arrangement. The oscillations of gi(r) is non-symmetric about the zero axis because of 
the typical columnar anti-ferroelectric structure of the simple cubic dipolar system. In a 
perfectly ordered state out of the six nearest-neighbours (r = 1) of a dipole two dipoles 
are parallel while four are antiparallel and therefore gi(r) is —1/3. For the case of next 
nearest-neighbours (r = 2) all six dipoles are parallel and therefore gi(r) is 1. For this case 
the long range behaviour of g2{r) (Figure [8b) is observed and its long range limit (0.427) is 
also found to agree with the square of the corresponding order-parameter (0.642) according 
to Eq. (jig). 

In Figure [8t we see that in the N v + phase (T* = 1.0) in spite of the presence of dipolar 
interaction (e* = 0.8) no anti-ferroelectric structure results. This is observed from the 
absence of long-range behaviour of g\{r). gi{r) remains zero for all values of r except within 
the first two neighbours where a small oscillation occurs. Also from the variation of g 2 {T) 
it is found that no long-range correlation of the transverse molecular axis is present which 
confirms the absence of biaxial order. Due to finite size effect a non-zero value (~ 0.1) of 
g2(r) is observed at large r which does not correspond to any true long range order for this 
case. 

Since in our study the dipolar part of the total pair potential depends on the orientation 
of the intermolecular vector, we have found, as expected from a previous study [3H], that 
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Figure 8: The first-rank and second-rank orientational correlation functions gi(r) (circles) 
and #2(7) (triangles) for the system of biaxial molecules having biaxiality A = 0.3 with and 
without transverse dipoles at the temperatures (a)T* = 0.3 (N B ), (b)T* = 0.8 (N B ) and 
(c)T* = 1.0 (Njj+). The dashed lines are used (in (6)) as a guide for the eye. 
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the directors are pinned along the lattice axes. The interaction energy of a pair of dipoles 
depends on the angle between the dipoles and the angles made by the dipoles with the 
vector joining them. The lowest energy for an interacting pair of dipoles occurs when they 
are collinear and directed along the vector joining them. The next minimum occurs when 
they are antiparallel to each other but perpendicular to the inter-dipolar vector. Therefore, 
for a simple cubic lattice system the ground state is a columnar anti-ferroelectric where the 
dipoles of a given row are parallel with one of the lattice axes and nearest-neighbour columns 
are antiparallel. Hence, due to the presence of anisotropy of dipolar interaction in orientation 
space the director frame coincides with the laboratory frame within a permutation of axes, or 
in other words, the principal and the secondary directors are pinned along the lattice axes. In 
order to verify the above fact we have evaluated the averaged quantities (B^ Q )i a b an d (R^iab 
which represent the uniaxial and biaxial order parameters with respect to the laboratory 
frame. The above quantities are calculated simply by taking the averages of Rq and i?22 
over the total number particles (i.e. the sample averages) for a given configuration and then 
by taking averages of these sample averages over large number of equilibrium configurations. 
The order parameters evaluated in the laboratory frame and their counterparts evaluated 
in the director frame differ by less than 2%. We have determined (R^^jiab for A = 0.3 and 
0.40825 and dipolar strength e* = 1.6 for different values of the temperature in ordered 
state. We have noticed that the directors remain pinned for all the temperatures and even 
close to the Nb <h- I transition they remain firmly pinned. 

5 Conclusion 

The phase diagrams of biaxial molecules possessing D 2 h symmetry and with transverse 
dipoles, along their shortest dimensions, have been studied using Monte Carlo simulation. 
As already elaborated our results are in good qualitative agreement with those in references 
[T5] and [2TJ. We find that the phase diagram changes significantly due to the influence 
of strong (transverse) dipolar interaction. The most significant effect of the presence of 
transverse dipoles is the splitting of the Landau point into a Landau line. Therefore, a 
range of geometrical structures of biaxial molecules are able to satisfy the most optimum 
condition of molecular biaxiality. We find that there is an optimal relative dipolar strength 
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for jj, ~ 3. 3D for which the width of the Landau line becomes maximum and for this optimal 
value the dipolar interactions have almost an equal contribution of energy with that from 
the dispersion interactions at low temperature. With further increase in dipolar strength the 
line begins to shorten and get shifted towards higher molecular biaxiality parameters. As 
the strength of the dipole is increased, the direct Nb -h- I transition temperature increases 
thus leading to an increase in the stability of the thermotropic Nb phase. The continuous 
nature of the direct Nb -h- I transition at the Landau point for the apolar biaxial system 
also appears to change because of the dipolar interaction. The sharp rise in the height 
of the specific heat peak and the sudden jumps in the order parameters at the Landau 
point due to strong dipolar interactions tend to indicate qualitatively that the nature of 
the transition becomes first order. A finite size scaling analysis is necessary to confirm this 
change in character of the phase transition. But considering the long-range nature of the 
dipolar interaction in simulations it becomes a very time consuming task and is therefore 
left out for future study. The structural properties in the uniaxial and biaxial phases are 
investigated by evaluating first rank and second rank orientational correlation functions. 
The dipole induced long range order of the anti-ferroelectric structure in the Nb phase 
is observed. The lattice model that we have used is very simple in comparison to real 
thermotropic liquid crystalline systems because real mesogenic molecules have flexibility 
[ID] and also possess both orientational and translational degrees of freedom. 

Finally we make a comment on the use of an anisotropic interaction (here the dipolar 
one) in a lattice model without taking an average over the orientations of the intermolecular 
vector. This is likely to result in orientationally ordered phases where the director becomes 
pinned and the structure of the phases depend on the lattice symmetry. Our observation of 
the appearance of the dipole induced long-range anti-ferroelectric order may have resulted 
from this. But we are not inclined to think that the biaxial nematic order we have observed 
in our work resulted from the use of a lattice model and the manner we have treated the 
dipolar interaction since the phases were identified by the relevant order parameters. A 
more elaborate study involving translational degrees of freedom of biaxial molecules with 
transverse dipole moment may perhaps be helpful. 
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